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Abstract 

We use the functional integral technique of Edwards and Lenard [pj to solve 
the statistical mechanics of a one dimensional Coulomb gas with boundary inter- 
actions leading to surface charging. The theory examined is a one dimensional 
model for a soap film. Finite size effects and the phenomenon of charge regu- 
lation are studied. We also discuss the pressure of disjunction for such a film. 
Even in the absence of boundary potentials we find that the presence of a surface 
affects the physics in finite systems. In general we find that in the presence of 
a boundary potential the long distance disjoining pressure is positive but may 
become negative at closer interplane separations. This is in accordance with the 
attractive forces seen at close separations in colloidal and soap film experiments 
and with three dimensional calculations beyond mean field. Finally our exact re- 
sults are compared with the predictions of the corresponding Poisson-Boltzmann 
theory which is often used in the context of colloidal and thin liquid film systems. 
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1 Introduction 



Up until 1961 the statistical mechanics of the classical one dimensional Coulomb gas 
was an unsolved problem. At more or less the same time the problem was solved by 
Lenard || and independently by Prager ||. A powerful alternative method of solution 
using functional integration was subsequently expounded by Lenard and Edwards |J. 
A good review of this work may be found in . It should be mentioned here that the 
two dimensional Coulomb gas may also be solved exactly at the temperature where 
e 2 /kT = 2; this exactly soluble case has been investigated in the electric double layer 
geometry by Cornu and Jancovici f7j. The problem of electrostatic interactions is one 
of profound importance in the theory of colloidal stability and also in the understand- 
ing of thin liquid films. In these problems one considers the behavior of an electrolytic 
fluid between two surfaces which either model the surface of large colloidal particles or 
the surfaces of the thin liquid film. The charging mechanism of the surfaces is usually 
of a statistical mechanical origin. For example in soap films made from sodium do- 
decyl sulphate (SDS), the soap anions have hydrocarbon tails which are hydrophobic 
and hence have a preference to lie on the surface of the film [p!2| . In colloidal sys- 
tems chemical reactions may occur between the colloid particles and the surrounding 
electrolytic medium again leading to surface charging. As the two planes are brought 
together surface charge regulation occurs. The precise qualitative behavior is still only 
understood within the context of mean field Poisson-Boltzmann type theories ||, ||, 
and at a more sophisticated level using the hyper- netted chain approximation (HNC). 
Rather surprisingly the HNC theory predicts, in the context of colloidal systems, that 
the electrostatic interactions between the planes may become attractive for small sep- 
arations 0, H]; this is supported by calculations of the fluctuations about the mean 



field solutions ||14|| . In the mean field model applied to soap films charge regulation is 
predicted |T2" |, but no attractive component is seen to appear within the electrostatic 
interactions. Interestingly the point at which charge regulation becomes important in 
the mean field model for SDS soap films coincides with the range at which collapse 



to a Newton Black Film (NBF) occurs [12] ■ There is much indirect evidence that the 



transition from a normal film to a NBF is of first order, for example it is believed to 
be exothermic and occurs via a nucleation process where regions of black film expand 
over the surface of the film. In this paper we propose to analyze the exactly soluble 



one dimensional version of the model proposed in [T2|]. We shall use the method of 
to solve the problem but we shall highlight the finite size effects appearing in the 
problem to gain an understanding of how charge regulation occurs in the model. We 
shall compare our exact results with those of mean field theory to ascertain, at least in 
one dimension, the accuracy of the traditional Poisson-Boltzmann mean field approach. 
The paper is arranged as follows. We formulate a form of the soap film model used in 
[ 12" ] in one dimension. The problem is solved using the functional integral formalism 
of (l| and the limit of bulk systems is rederived for the sake of completeness. We 
then analyze the problem in the case of finite films with surface binding interactions 
and discuss the nature of the charge regulation and the stability criterion for the one 
dimensional film. We then compare the mean field Poisson-Boltzmann theory with 
the exact results. Finally we conclude with a brief comparison between the qualitative 
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behavior observed in the one dimensional system and that of experiments. 



2 Analysis 

Here we shall summarize the approach of |lj and apply it to the system in which we are 
interested. The field theory for the system we shall consider is derived from considering 
a model consisting of a monovalent soap molecules whose anions are attracted to the 
surface of the soap film by the presence of an effective potential V(x) which acts on them 
and whose support is localized at the two adjacent surfaces of the film. In addition one 
may add an additional monovalent electrolytic species. In the grand canonical ensemble 
if the fugacities for the soap anions/cations and the electrolyte anions/cations are given 
by n s and \x e respectively, then the partition function is given by [l] 




d[if)] exp 



(d^ 
\ dx 



dx 



+ /i s I (exp (— i/3 2 ei/j — f3V(x)J +exp (i^eip^ dx + 2^i e I cos (j3*ei/)j dx (1) 

The derivation of the representation above is quite standard and relies on the intro- 
duction of the Hubbard Stratonovich field ip. In general dimension the field theory is 
a form of the Sine-Gordon field theory and is not generally soluble. For simplicity we 
shall chose the form of V(x) to be such that 

exp(-f3V(x)) = 1 + X*8(x) + X*5(x - L), (2) 

i.e. the effective surface potential is highly localized about the boundary points and 
L. The A 
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in 



appearing in (0) is similar to the adhesivity introduced by Davies 
his analysis of the surface tensions of hydrocarbon solution although the idea of such 
a surface active term goes back to Boltzmann. For physically realisable soap films 
V(x) is not strictly localized as the effective potential created due to the hydrophobic 
nature of the soap anion hydrocarbon tails has a support over a region of the length 
the tail between the surface and the interior of the film (see IHH for a discussion of 



the mechanism generating this potential). However for the purposes of demonstrating 
the the essential physics of charge regulation in a one dimensional system our choice 
of V(x) should be adequate. With this choice of V we obtain 





— — dx 



d[ip] exp 



+ 2/i / cos {j3^eip^j dx + A (exp(— i(3* eip{tij) + exp(— i(3*eip{L)) 



(3) 



with jj, = fi s + n e and A = /z s A*. Because now the potential V acts only on the end 
points we may write Z in path integral notation as, 

2 



1 

N 

+ 2/i 



dipod^L I d[ip] exp 

V'o 




dx 



cos 



(/3*eV>) dx + A (exp(-i(3^eip(0)) + exp(-i/3^eip(L))^ . (4) 
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The path integral is that for a diffusing particle in a cosine potential (the normalization 
coming from the determinant is absorbed in the factor N), consequently following [l[ 
we find 

1 f 11 
Z = — dijj dip L P (ip ,ip L \L) exp A (exp(-i^ ei/) Q ) + exp(-i/52e^L)) , (5) 



N 

where P(i[)o,i/;\x) obeys 

1 d^P dP 

__ + 2MC08 (^)P=-, (6) 

subject to the initial condition P(ip ,ip\0) = 5(ipo — 4>)- Note here that the boundary 
terms ipo and ipL in our path integral are free and are integrated over, this is in contrast 
to the study of [Q, where ipo — and i/jl is left free. This is because in the formulation 
of 0, it is assumed that the overall system bulk plus the film is electroneutral and that 
the bulk lies to the left of the point 0. In our case the bulk may be screened from the 
film and the most general boundary conditions are those we have employed. In order 
to determine the normalization factor iV we note that when e = A = then we should 
obtain the ideal gas result Z = exp(2/iL). In this case 

P(ip , ip\x) = 1 , exp (-— (tp - tp ) 2 ) exp(2/ix) (7) 

giving simply 



Z\ e=x =o = — J <#oexp(2/iL). 



At this point we must regularize the partition function by bounding the possible values 
of ipo between two extrema yielding 

N= (l ^ ■ (,Y ^ 

max(^oj — mm^oj 

In the case where e ^ we may use the fact that the action is invariant under transla- 
tions of 27r/e/32 to obtain (we have simply assumed that the extremal values of ipo are 
integer multiples of 2-n/e(3^) 



Z = — — / dtpo / d?p L exp A (exp(-i/32e^o) + exp(-i/3^e^ L )) P (*jj ,ip L \L) . 
Zir Jo J-qo 1 v /J 

(10) 

A further simplification is obtained by noting that 

P (Vb, ML) MlWl = K(ip Q \L), (11) 



where K(ip \L) obeys 
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with 



# = ^1T72 + 2 A tcos (/^o), (13) 

subject to the initial condition K(ip \0) = fijpo)- In our problem 

f(ipo) = exp (A exp(-i/35e-0 o )) (14) 

For simplicity in notation we shall take (3 = 1 from here on. To recover the temperature 
dependence the rescalings e 2 — > f3e 2 and P — > /cTP should be performed. The final 
expression for Z is 



2?r Jo 

which in operator notation may be expressed as 



/ e # i^ |L)/(^o), (15) 
Jo 



#o/(^b) exp(Lif)/(^). (16) 



2n Jo 

The free energy of the film is simply given by .F(.L) = — log(Z). In order to calculate 
the disjoining pressure we assume that the film is attached to an infinite reservoir of 
particles with the same chemical potentials. In addition the over all volume (i.e. length) 
of the system is conserved. This last point is quite important in understanding the 
statistical physics of small systems where thermodynamic ideas cannot be necessarily 
applied directly. Our physical assumption of the over all incompressibility of the system 
is motivated by the principle that the background space of water can be assumed to 
be incompressible in most experimental regimes in colloid and thin film science. 
Therefore on changing the volume of the film by 5L, the free energy of the bulk is 
changed by —5L0OO where doo is the free energy per unit volume of an infinite system 
(i.e. one which is really in the thermodynamic limit) without any boundary interaction. 
The work done on this change is PdSL where P<i is the disjoining pressure. Therefore 
we find that 

dF 

Pd = ~dL + 6 °°- (17) 
Finally defining the field = efJ^ip we obtain 

Z=±- I*" dfofifo) exp(L#)/(0 o ), (18) 

271 JO 

where now f((f>o) = exp(Aexp(— i<f>o)) and 

^ = e 2 ^ + 2 /U cos(0 o ). (19) 
One point observables are calculated as 

1 r 2 ^ 

(0(<P(x))) = -J Q # o /(0o) exp(xH )O(0 O ) exp((L - x)if )/(<fo), (20) 
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the extension to n-point functions being trivial. Here for convenience and to highlight 
the different regimes we develop the notation. We write 

e 2 

H = —H* (21) 
2 v ; 

where 

H* = 1 -^ + acos( ( p ) (22) 



with a = 4/i/e 2 . 



3 Results for Large Films: Thermodynamic Limit 

The eigenfunctions of H periodic on [0, 2ti] are the periodic Mathieu functions Xn((fi, a) 
whose eigenvalues we denote by 7 n («) and where it is easy to see that the largest 
eigenvalue 70(a) < a. Hence in the case where e 2 L/2 <C 1 

Z ~ h ( J ( 2 ^^o(0)/(0)) 2 exp(^L 7 o(a)). (23) 

therefore = — y7o(a) and any boundary terms become insignificant in the thermo- 
dynamic limit. The bulk pressure is then Pbulk = —Goo ■ 



3.1 Small a , large e 

In the case where a is small one may evaluate 70(a) in perturbation theory and we find 
that M 



1 2 7 , 29 fi 68687 R 

7o (a) = -a 2 a 4 H a 6 a 8 

' v ; 2 32 144 294912 

In this regime the average density is extensive and is given by 



+ 



P 



dP 



e d , yo(a) 

— a — ; — 

2 da 



The condition that a is small is therefore equivalent to 

^2~p/e « 1 . 



a ~ 



(24) 



(25) 



(26) 



This implies that the electrostatic energy is much greater than the contribution from 
the entropy. From equations (|17 , |24|, ^5|) the pressure is then given by the series 



bulk 



P 



1 + 



7 



23 
288 



4897 fp 
122288 Ve 2 



+ 



(27) 



This result was explained by Lenard as an effect of dimerization. The leading term is 
independent of e and is the perfect gas result for a density of p/2. This is explained by 
the positive and negative charges binding in pairs to give, in leading order, a neutral gas 
with half the original particle density. The non-leading terms correspond to multipole 
interactions such as van der Waals forces etc. 
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3.2 Large a , small e 



The calculation of 70(a) for large a can be formulated as a perturbation series for H* 
in equation (|22"D obtained by expanding the cosine and writing 



H* 



a + H osc + 0(a< 



:28) 



where H osc is the harmonic oscillator Hamiltonian 



d 2 



-a (f) 2 



(29) 



Thus the perturbation theory is for the anharmonic terms in equation (^) using the 
basis of oscillator states associated with H osc . The first term in the pressure is due 
to the 0(a) term in equation ( |2"8| ) and gives the free gas contribution. The next 
correction arises from the ground state eigenvalue of the oscillator and is 0(\fa) . This 
is the well-known Debye-Hiickel term. We expect a power series in a _1//2 , but to carry 
out the perturbative expansion becomes increasingly difficult as the order increases. 
Instead, we can formulate the problem as an expansion in Feynman diagrams. A similar 
approach to the calculation of the electrostatic free energy of a system of fixed charge 
macroions has been used by Coalson and Duncan fl8j and Ben-Tai and Coalson |19| . 
In the bulk limit we write the Feynman kernel for as 



JC(L, a) = N J dip exp (j dx 



1 / dip\ a 
"2 U + 2 C ° SW 



(30) 



where iV is the normalization factor chosen so that /C(L, a = 0) = 1 . In the limit 

L — ► 00 

/C(L,a) = exp(^L) => 7o (a) = 2^ log/C(L, a) . (31) 

It is convenient to discretize L so that L = ne , where e is the lattice spacing. Then 
the operation d/dL can be performed directly on the LHS of equation (RO) as 



1 d 

70(a) = 2-— log /C(L, a) 
n oe 

= «y) 2 > + a(cos(V0> 



(32) 



The last term arises because the normalization factor iV depends on e . This term 
cancels trivially with a simple e divergence in the 1-loop graph for ({ip') 2 ) ■ We can 
now take the limit e — > . The density, defined by equation (p5|), is given by 



C 2 CL 

P = -y(cos(V>)) 



(33) 



and thus the second term in equation ([32] ) corresponds to the free gas term. Hence, we 
have 



bulk 



e 



m 2 ) 



(34) 
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We define the Debye mass m by m 2 = a/2 and then 

K{L,a) — J # exp(-S(i/))) , (35) 
where overall irrelevant constant factors have been omitted, and 

"1,,a2 1 2,2 (-l) n+1 m 2 



S(ip) = dx 



oo 



oo 



(iP'Y + » 2 + v Z \T~~ ^ 



(36) 



2 vr y 2 r (2n)! 

A standard Feynman graph expansion of closed loops for {{ip') 2 ) can be obtained and 
hence the pressure can be calculated from equation (34). Standard dimensional analysis 
shows that the series obtained is in inverse powers of m and that a diagram with I loops 
behaves as m 2 ~ l . To three-loop order we evaluate the diagrams shown in figure [I] and 
find 

e 2 / m 1 3 \ 

P - = " + 2 (" 2 + 16 + 512^) (37) 
In order to express Pbulk as a function of the density, p, we use equation ( |33"D and 
calculate p in the loop expansion. As before, the diagrams with / loops behave as 
m 2 ~ l . If P is calculated to Z-loop order, then p is needed to (/ — l)-loop order. To 
two-loop order we find 

" - eV i 1 ~ sb) ■ w 

Note that, alternatively, 

p = y^T' (39) 

and if we define t = log(m) and write Pbulk = P + P e then we have 

dp dP e , . 

which has solution 

p r°° 1 dP e , . , 

—2 = / — 2 dm + 1 , 41 

Jm m z dm 

where we have used the boundary condition that p/m 2 — > 1 as m — > oo . 



Then P(m) can be re-expressed as a series in yp/e 2 : 

1 / 1 /e^ 

Pbulk = P - -\/pe 2 + \ — + ■■■ ■ (42) 

p 4 VF 1024 V p 

This agrees with Lenard and it is relatively easy to evaluate the loop expansion to 
higher orders to improve on Lenard's result. The second term is the familiar Debye- 
Hiickel contribution and it should also be noted that the two-loop contribution is zero. 



3.3 The bulk pressure 

In figure ^| we show the computed value for Pbulk compared with the predictions of the 
previous two sections. For convenience Pbulk has been scaled by a factor of e 2 /2. As 
can be seen the curves from equations ([24] and ( 37|) fit very well except in the region 
0.7 < a < 1.1 where even so the discrepancy is not very large. 
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4 Thin Films: Finite size Effects and Surface Charge 
Regulation 

When the intersurface distance of the film becomes small in the sense that \e 2 L is no 
longer large, then we may not apply the thermodynamic result (|23|). 
However if |ae 2 L and \e 2 L are both small, which is certainly the case for extremely 
small L, then one may expand the operator exp(LH) in powers of LH. To second 
order in LH one obtains 

L 2 

Z = 1 + 2fiLX + — (2 / u 2 (A 2 + 1) - /iAe 2 ) + 0(L 3 ) (43) 

thus yielding F ~ — 2/iLA and Pa = 2fiX + 6^. Hence one has the limiting value of the 
disjoining pressure is negative if the value of A is sufficiently small. However, one has 
the bound that 70(a) < a and hence 6*00 > — 2/i thus 

limP d > 2//A-2//. (44) 

Hence the film certainly has a positive disjoining pressure at small differences if A > 1. 
The stability of the film at small separations is determined by 

BP 

lim ^ = 2// 2 (l - A 2 ) - /iAe 2 . (45) 

If this term is positive then the film collapses to the point thickness L = 0. For this to 
happen one must have 

A < 1, (46) 

and 

" > 2(1^ < 47 > 
The value of the surface charge a is given by 

a = -Ae(exp(-i0(O))) = -eLXfi . (48) 

Hence over short distances the surface charge decays linearly as the two surfaces are 
brought together. 



5 Intermediate Regime 

In the regime between very thick and very thin films we shall resort to a numerical 
analysis of the problem. There are two methods of interest which we detail in the 
following sections. For convenience of notation we shall work in units where L is scaled 
by e 2 /2. That is, a factor of e 2 /2 is absorbed into all length variables. 
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5.1 The Mathieu function method 



The disjoining pressure and other properties of the film can be calculated using the 
even and odd Mathieu functions that are the eigenf unctions of H* defined in equation 
(p2j). The kernel /C(L, a) defined in equation (|30| ) can be computed as an expansion on 
the Mathieu function by resolution of the identity on the basis of these states. In this 
way the disjoining pressure may be in general be written as 

p m _ e 2 E^Li (7o («) ~ 7n(a)) c 2 n exp(7 n (q)L) 
d[ } " 2 £~ cl exp( 7n (a)L) ' [ ^ 



where 

r 2n 







d0x n (0)exp(Aexp(-«0)). (50) 



If the eigenvalues of H are arranged in descending order i.e. 70 > 7i > 7 2---> the 
corresponding eigenfunctions are even about it for n even and odd about tt for n odd. 
Hence c n is purely real for n even and purely imaginary for n odd. If A = then 
c n = for n odd and hence Pd is always negative, hence the force between the two 
interfaces is always attractive. One sees from the above expression that it is the even 
wave functions which are attractive and the odd wave functions which are repulsive (as 
the demoninator on the RHS of ([EJ) is Z and hence positive). At long distances 

2 2 

p d ~ %r\ (7o - 7i) exp(L( 7l - 7o )) > (51) 
/ c 

Hence for a non-zero A the long distance disjoining pressure is always positive since 
c\ < . It is clear however that the disjoining pressure may become negative at smaller 
values of L. 

The anion and cation number densities as a function of x, the distance through the 
film, may also be calculated. Denoting these densities respectively by p + and p_ we 
find 

P± = Mf*\Xn) exp(-7 n x) (Xn\e ±l ^\x m ) exp(~7 m (L - x)) (Xm\f) , (52) 

where 

(Vlx) = / d<W0M0) • (53) 



We are able to construct both the even and the odd Mathieu functions and their 
eigenvalues for any value of a using Given's method for diagonalizing a tri-diagonal 
matrix. The eigenfunctions of H* are found on a discretization of the interval [0, 27r] and 
the appropriate matrix elements in equations (f49|, |52"|) can be calculated numerically. 



5.2 Fourier method 

It turns out that there is a more direct method to calculate the disjoining pressure which 
exploits the periodicity inherent in the system. This method is especially effective for 
low temperature (small e). It does, however, become much more complicated when 
other observables such as the density profiles are being calculated. Expanding K in 
terms of its Fourier modes, i.e. writing 
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oo 



K((f)\x) = exp(ax) 2J b n (x) exp(in0), (54) 



n=— oo 



one has that 



A _n 

6«(0) = t rrforn^O 

— n ! 



for n > (55) 



and the b n evolve via the equation 

db 



= - n 2 b n + a(b n+1 + 6 n _! - 26 n )/2. (56) 

dx 

Finally the partition function is given by 

oo \n 

Z = exp(aL) £ -— 6 n . (57) 

n=0 W' 

In this regime we shall also be interested in the mean value of the surface charge a. 
Including a source term in the original formulation of the problem it is a simple matter 
to show that 

Ae 1 r 2n 

a = Ae(exp(-z0)) = — — / dcj) exp(-i(j) )f(<p )exp(LH)f((f) ). (58) 

Z Ztt Jo 

In terms of the Fourier expansion this becomes 

■y oo 

a = — exp(aL) ^ --^b„ n+1 (59) 
L n=o v 1 )- 

The disjoining pressure may be computed similarly. 

In what follows we shall consider three cases which are paradigms for the different 
regimes of high, intermediate and low temperature. Since we have set kT = 1, high 
T corresponds to a small charge parameter, e and vice-versa. Apart from an overall 
dimension-carrying factor the results depend on e and /i through the combination 
a = 4fikT/e 2 , and in what follows we choose // = 1 and hence in our units a = 1/e 2 . 
The three regimes of temperature are characterized by the three values of charge: 
e = 0.1,1,4 . 



5.3 e = 0.1, a = 400 

From equation fl5?D the bulk pressure is P = 1.926. The major correction to the free 
particle pressure, P/ ree = 2, is the Debye-Huckel term and the two and three- loop 
contributions are a correction of only AP = 0.003 . In figures ||| and § we show the 
pressure P versus film thickness L for various values of A in the range 0.9 to 1.2 . Also 
plotted is the prediction for the bulk pressure to which all curves should be asymptotic. 
As can be seen there is a collapse in all cases shown for A. The details of the collapse 
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differ, however, as A increases. For the lower values of A the collapse is to a film of 
zero thickness which would, of course, be dominated by the detailed structure of the 
surface physics which we have subsumed in to a layer of zero thickness. Two maxima 
are clearly visible for A = 0.93, 0.95. The one at larger L is the location of the ordinary 
collapse point. The maximum at smaller L and the consequent multiple- valuedness of 
the curve in L versus P in this region implies a hysteresis phenomenon as P is cycled 
for very thin films. This kind of effect is reminiscent of a first-order transition which 
predicts that for a 3D film there will be domains of different thicknesses which will 
grow or contract like 2D bubbles. Of course, it remains to be seen whether intuition 
from ID survives for the realistic 3D case. The typical length scale is 0(e 2 /2) ~ 5 10~ 3 , 
which is very small compared with the values of L plotted in figures fj and (|. 
For the larger values of A plotted the curves the collapse is to a thinner film but not 
to one of zero thickness, as A increases the maximum at small L eventually disappears 
and for much larger A the collapse phenomenon itself disappears. 
In figure we shown the surface charge a defined in equation (|58|) as a function of L. 
There is no feature which hints at the presence of the collapse phenomenon appearing 
in the associated pressure curves, but in all cases a decreases with L. For small L the 
behavior agrees well with the prediction of equation (f48[). 

The anion and cation densities have been computed as a function of x for various 
values of L using equation (|52|) . For e = 0.1 the variation with x is mild and shows no 
features of note. We show the midplane values for each species as a function of L and 
for various values of A in figure ||. 

It is interesting to note that both methods described in subsections |5.1| and |5.2| were 
used to calculate the disjoining pressure. However, while 10 fourier modes were amply 
sufficient, the number of Mathieu modes needed was 40. In particular this large number 
of modes was found necessary to reproduce the secondary collapse maxima shown in 
figure H 

5.4 e = 1.0, a = 4 

As in the previous section the pressure P versus L is plotted in figure ^ for values 
of A in 0.3 < A < 0.8 which span a region of collapse. In this case there is just 
one point of collapse to a film of zero thickness (in our approximation) and which 
disappears for A between 0.7 and 0.8. From equation fl3"T| ) the bulk pressure is predicted 
to be Pbuik — 1-3262. From the exact calculation we find Pbuik — 1.32584 which is 
in good agreement with the prediction. To guide the eye the computed asymptotic 
value is shown in figure [7j. These pressure curves are well reproduced by the Mathieu 
function method with as few as 8 modes. Unlike the case in the previous section 
a = 4 is sufficiently small that the physics is dominated by the lowest-lying Mathieu 
eigenf unctions. This is mainly due to the fact that A is smaller and so the overlap 
(Xn\f) falls °ff more sharply with n. This in turn means that the pressure peak occurs 
for larger L than in the a = 400 case. However, the large-a result for the bulk pressure, 
equation ([IT], , still holds very well in this region which means that the Debye-Hiickel 
approximation is good. 

The surface charge a is plotted in figure |8] and, as in the previous case, there are no 
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features associated with the pressure maxima of figure [7] 

The anion and cation number densities as a function of distance, x, through the film 
are shown in figure ^| for A = 0.5 . These quantities were calculated using equation 
(|52]). The anion (cation) curves are the higher (lower) set in this figure. There are no 
unusual features and the curves for the other values of A in 0.3 — > 0.8 are of similar 
form. 



5.5 e = 4.0, a = 0.25 

The pressure P is plotted versus L for A = 0.12,0.121,0.122,0.123 in figure |T0|. The 
collapse region is again evident but it should be noted that it occurs only for a very 
narrow range of A values. Of course, A is a parameter that is determined by other 
variables and is not fixed externally. The bulk pressure is no longer given by the large- 
a expression (equations ( |37| |4~2|) but is well fitted by the small-a result (equations 



27]). The prediction is Pbuik = 0.243529 and the computed value is Pbuik = 0.243531. 
This value is shown for reference in figure |TTj| . 

The curves for the surface charge, a, are similar to those of previous sections and are 
not reproduced here. The small L behavior is again consistent with equation (|41j|). 
The anion and cation number densities (equation ( |52]) are plotted for A = .123 and 



L = 0.1,0.51.1 in figure [11]. It is interesting to note that for both species the density 
falls sharply at the film surface and the anion density reaches a peak for the thicker 
films which is located only a short distance into the film. The position and shape of 
this peak is independent of thickness L and seems to be a universal feature of the low 
temperature case. Only the lowest 4 Mathieu modes make an appreciable contribution 
since the effective charge is large and from equation (|52] ) this causes a strong exponential 
suppression on all but the lowest modes (note that in equation ([52]) a factor of e 2 /2 is 
absorbed into all lengths). Also, the values of A in the collapse region decrease as e 
increases and so the surface function /(</>) (equation ([14]) ) oscillates less fast and only 
has appreciable overlap with the lowest modes. For these reasons the species number 
densities are dominated by the contributions from the lowest modes and so show more 
structure at low temperature than at high temperature. This is to be expected since the 
electrostatic energy dominates the thermal energy. We have observed similar maxima 
in the density profiles for other values of e if A is sufficiently small. This is due again 
to the dominance of only a very few low-lying Mathieu modes. 



6 Comparison with Poisson-Boltzmann Theory 

The Poisson-Boltzmann (PB) theory for our system may either be derived directly by 
standard thermodynamic techniques || ||, or as the mean field theory for the field 
theory ([I]). The theory has been used in a wide context in soft condensed matter 



physics and in particular to analyze the behavior of soap films in [jlT| , [12] , and also 
in the context of colloidal stability ||. In general it is fair to say that it has been 
reasonably successful in predicting the physics of systems where interplane distances 
are reasonably large and for monovalent ionic species |T5] |. 
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The resulting equations are (again scaling so that (5 = 1) 



r2 J 

-% = 2//esinh(e0) + Xe(5(x) exp(e0(O)) + 6(x - L)) exp(e0(L))), (60) 
dx 2 

where here is the mean field electrostatic field. Assuming symmetry about the 
point L/2 (however see the comments in the conclusion) and using the condition of 
electroneutrality (which is a mean field assumption), the boundary conditions are 



dx °' L dx 1 2 



L 



(61) 



Interestingly <p appears as a purely imaginary saddle point of the theory ([!]). In the 
region [0 + , L~] the above equation reduces to 

= 2/iesinh(e0), (62) 

with the boundary conditions 

%=0 (63) 

dx 2 

and 

^| o = -Aeexp(e0(O))=a, (64) 

where a is the surface charge. It easy to show that the mean field free energy over the 
bulk is H 

F MF = / ((f)') 2 dx - (2Aexp(/5e0(O)) + 2,aL(cosh(e0(L/2)) - 1)) . (65) 
J o 

Then we find 

P d = 2 /U (cosh(e0 m ) - 1) , (66) 

where m = 0(L/2) is the midplane potential. One immediately sees that in the case 
A = then = is a solution and the film is always marginally stable in the sense that 
Pd = 0. In general any non-zero A gives a non-zero value of m and hence the film is 
always stable for non-zero A. This is clearly at variance with the exact results derived 
here. Moreover, the mean field bulk pressure is Pbuik — 2/i which is only applicable to 
the limit e — > or, equivalently, a — > oo. 

In general one must resort to a numerical solution of the above mean field equations. 
However in the case where L is small such that <fi varies only slightly we may use the 
approximation, 

= m + C(a;-^) 2 . (67) 
Substituting this into equation fl62"D yields C = fie sinh(0 m ). Using this in the boundary 



condition (|64 ) then yields 

/iL sinh(e0 m ) = — Aexp(e0 m ) + 0(L 2 ). (6£ 
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Solving this yields 



1 9 A 

-^M^ + l). (69) 



Hence in this limit 



and 



P d (L)~\ 2 -^ (70) 



One sees that while the surface charge a does decay to zero it does so as in com- 
parison with the exact result (EST) . In addition the disjoining pressure ( [70]) actually 
diverges rather than tending to a constant. 

For infinitely thick films one may use the condition that <p m — > in order to calculate 
the surface potential. In this case the surface potential is given by the equation 

cosh(e0(O)) - ^-A 2 e 2 exp(e0(O)) = 1, (72) 
from this we find that the physical solution is 

0(O) = -^log(l + Ae/^) (73) 

giving a surface charge 

- = w ^ (74) 

1 + Ae/v^Ti v ' 



At intermediate distances one has to numerically solve the PB equations. For the cases 
discussed in section |5| there is no agreement at all between the numerical solution to 
the PB equation and the exact result. This is to be expected since there is no collapse 
predicted by the PB equation. However, there is no agreement even on the rising part 
of the pressure curve at L much greater than that at the pressure maximum. Also, the 
values of Pbuik are not close to the mean-field prediction of P^K = 2-0 although for 
e = 0.1 this value is approached. Nevertheless, in this latter case there is still a large 
disagreement between the exact and mean-field curves. Indeed, we have investigated 
very small values of e for a large range of A values but have not found any reliable 
agreement between the exact theory and the PB equation. The PB equation may 
be applicable for even smaller values of e than we have investigated. Indeed, a naive 
analysis of the applicability of the saddle point method for the theory (|1|) suggests that 
a should be large, i.e., either fi ^> 1 or e 2 <C 1, thus giving either \i or 1/e 2 as the large 
parameter justifying the saddle point analysis. However, in the cases we have analyzed 
here, mean-field theory and the PB equation are of very little value in the analysis of 
the one-dimensional Coulomb gas. 
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7 Conclusions 



In conclusion we have derived an exact solution for the one dimensional Coulomb gas 
with boundary effects. Surprisingly the mere presence of a boundary, without any 
surface adhesion term, leads to a reduction of the density near the boundary with 
respect to the bulk. This effect means that the disjoining pressure of the system is 
negative and the resulting film will tend to collapse. When A > we have shown that 
at sufficiently large distances the disjoining pressure must be positive, and hence a 
stable common film regime exists. However, if the value of A is not too large a collapse 
phenomena may occur where the disjoining pressure decreases as the surfaces come 
together. The disjoining pressure may even become negative signalling the onset of 
strong attractive forces in the system; this may well be the one dimensional version of 
the collapse to a NBF seen in experimental systems. We have also seen the possibility 
of secondary collapses in the parameter ranges we have studied; it would be interesting 
if one could find an experimental system exhibiting a secondary collapse. In principle 
multiple collapses are possible, but we have yet to see more than two. 
Poisson-Boltzmann theory predicts a stable film for any non-zero value of A and in 
addition the calculated mean field disjoining pressure is larger that that of the exact 
calculation. Taking into account the full theory and all its correlations does indeed 
introduce an attractive interaction over and above the mean field result, in accordance 
with the calculations made in three dimensional systems using techniques beyond mean 
field theory ^ |^, [TJ| . We would like to comment here that in our solution of the mean 



field equations we have, as is done throughout the literature, always assumed that the 
mean field solution is symmetric about the midplane of the film. In physical terms this 
seems quite plausible for thick films where the two planes do not interact and hence 
there can be no breaking of spatial symmetry. The variant of mean field theory used 



in [12|, uses this symmetric solution and the theory describes extremely well both 
surface tension data for SDS bulk solutions and the disjoining pressure isotherms up 
to the point where the collapse occurs. One may show [T^| in the field theoretic sense 



that the mean field solutions we have found here and in |T2|, [13] are indeed stable local 
minima of the free energy and hence they at least describe a metastable state. The fact 
that the mean field solutions work so well in this context up to the collapse point suggest 
that another mean field solution with a broken spatial symmetry and possibly with a 
complex part may appear with a lower free energy than that of the symmetric real 
solution. This would also be consistent with the experimental indications (and indeed 
conclusions that may be drawn from our exact solutions in 1 D) that the transition to a 
Newton Black Film is of first order. Work on this problem is currently under progress 
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3-loops 



Figure 1: The one- dimensional Feynman graphs up to three- loop order which con- 
tribute to the calculation of the bulk pressure in equation (f?3|). The operator insertion 
for (ip') 2 is shown and corresponds to the insertion of the factor p 2 in the appropriate 
loop integral. 
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Figure 2: The bulk pressure Pbuik, scaled by e 2 /2 versus a = 4/i/e 2 . The solid line is 
the exact computed curve and the dashed and long-dashed curves are the predictions 
of equations ( p4]) and fl37D , respectively. The predictions fit very well except in the 
region 0.7 < a < 1.1 . 
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Figure 3: The pressure P versus film thickness L for kT = 1.0, e = 0.1 and fi = 1.0. 
for < L < 10.0. The different curves are for A = 0.9, 0.93, 0.95, 0.97, 1.02, 1.2 which 
respectively correspond to the curves from lowest to highest pressure at any given L. 
The phenomenon of primary and, in some cases, secondary collapse are clearly visible. 
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Figure 4: The curves shown in figure y extending to larger L to highlight the region 
of primary collapse. As in figure [| the higher the value of A the higher the pressure at 
given L. 
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Figure 5: The surface charge a defined by equation (fl8l) versus film thickness 
L for kT = 1.0, e = 0.1 and /i = 1.0. The different curves are for A = 
0.9,0.93,0.95,0.97,1.02,1.2 corresponding to the curves from lowest to highest a. 
There is no feature which hints at the presence of the collapse phenomenon appearing 
in the pressure curves shown in figures [3] and |j. 
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Figure 6: The midplane anion (upper curves) and cation (lower curves) densities as a 
function of film thickness L, for e = 0.1, /x = 1, A = 0.9 (solid) and 1.2 (dashed). 
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Figure 7: The pressure P versus film thickness L for kT = 1.0, e = 1.0 and fi = 1.0. 
The different curves are for A = 0.3, 0.4, 0.5, 0.6, 0.7, 0.8 which respectively correspond 
to the curves from lowest to highest pressure at any given L. The collapse phenomenon 
occurs for the smaller A values and disappears between A = 0.7 and A = 0.8. 
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Figure 8: The surface charge a defined by equation (fi8"D versus film thickness 
L for kT = 1.0, e = 0.1 and \i = 1.0. The different curves are for A = 
0.9,0.93,0.95,0.97,1.02,1.2 corresponding to the curves from lowest to highest a. 
There is no feature which hints at the presence of the collapse phenomenon appearing 
in the pressure curves shown in figures [3] and |j. 
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Figure 9: The anion (solid curves) and cation (dashed curves) densities as a function 
of distance x through the film for e = 1.0, fx = 1, A = 0.95 and for various film thickness 
L. 
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Figure 10: The pressure P versus film thickness L for kT = 1.0, e = 4.0 and [i = 1.0. 
The different curves are for A = 0.12, 0.121, 0.122, 0.123 which respectively correspond 
to the curves from lowest to highest pressure at any given L. The collapse phenomenon 
occurs for the smaller A values and disappears between A = 0.122 and A = 0.123. 
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Figure 11: The anion (solid curves) and cation (dashed curves) densities as a function 
of distance x through the film for e = 4.0, /i = 1, A = 0.123 and for film thicknesses 
L = 0.1,0.5,1.1. 
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